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Abstract. The results from an X-ray spectroscopic study of the giant elliptical galaxy NGC 5044 in the center of a galaxy 
group are presented. The line dominated soft X-ray spectra (mainly Fe-L and O VIII Lya) from the diffuse gas are resolved for 
the first time in this system with the Reflection Grating Spectrometers on-hoard XMM-Newton and provide a strong constraint 
on the temperature structure. The spectra integrated over 2' 20 kpc) in full-width can he described by a two temperature 

plasma model of 0.7 keV and 1.1 keV. Most of the latter component is consistent with originating from off-center regions. 
Compared to the isobaric cooling flow prediction, the observation shows a clear cut-off below a temperature of 0.6 db 0.1 keV. 
Furthermore, the Fe and O abundances within the central 10-20 kpc in radius are accurately measured to be 0.55 ± 0.05 and 
0.25 ±0.1 times the solar ratios, respectively. The observed cut-off temperature of this galaxy and other central galaxies in 
galaxy groups and clusters are compared with the scale of the galaxy and properties of the surrounding intra-cluster medium. 

Based on this comparison, the origin of the lack of predicted cool emission is discussed. 
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1. Introduction 

At the center of X-ray luminous galaxy clusters and groups, 
a giant elliptical galaxy is often found. The formation of the 
galaxy most likely involves a series of galaxy mergers. Around 
these regions. X-ray emission is also peaked, indicating a con¬ 
centration of the hot gas (one to a few keV). The central gas 
density (> 10^ m“^) may result in gas cooling along with mass 


in-flow (cooling flow; CF; e.g. Fabian 1994), although the pre 


dieted flow velocity has not been directly measured, and the 
fate of cooled gas has not been established. Instead of in-flows, 
the cooling may be counterbalanced by heating mech anisms , 
such as thermal conduction (e.g. Takahara & Taka hara 1979 ), 
magnetic reconnection (e.g. Soker & Sarazin 
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detailed spectroscopy with A5CA, XMM-Newton, and Chandra 
in representative clusters revealed a lack of soft X-ray emis¬ 
sion as predicted by a simple CF picture (e.g. Ikebe et al. 


1999; Peterson et al. 2001; David et al. 2001; but also see e.g. 


199C; Makishima 


et al. |2001[ ), or AGN activity (e.g. Churazov et al. [2001[ ). In any 
case, this cooling, mass in-flow, and heating processes should 
involve energy and/or mass transfers on the largest scale in the 
universe. 

X-ray observations established the presence of plasma near 
the cluster center considerably cooler t han th at found further 
out in the cluster (e.g., Canizares et al. 1979 ). However, a se¬ 
ries of innovative X-ray results on nearby clusters have cast 
a serious doubt on the reality of CF hypothesis. In particular. 


Schmidt et al. 2001 ). The detailed spatial and thermal structure 
of this cool component has not yet been systematically mea¬ 
sured. Consequently, the origin of the lack of cool emission is 
not fully understood. 

Here we report the RGS spectroscopic observations of the 
core of the NGC 5044 galaxy group (WP23). X-ray observa¬ 
tions showed that the X-ray emitting gas is distributed over the 
group and c entere d on the giant e lliptic al galaxy, NGC 5044 
(David et al. 1994; Fukazawa et al. 199f). The gas temperature 
is about 1 keV with a central drop down to 0.8 keV. Compared 
to these ROSAT and ASCA observations, the RGS provides 
much better spectral resolution in the 0.3-2.0 keV band; at 
1 keV these instruments have energy resolution (E/AE) of 
~ 2,^ 10, and ~ 200, respectively. This merit of the RGS 
overwhelms its two disadvantages for extended sources, that 
it cannot yield two-dimensional distributions of the emission, 
and that the source extent degrades the energy resolution. We 
present here one of the best quality soft X-ray spectra of el¬ 
liptical galaxies, and use the spectra to constrain the thermal 
structure and metallicity within the central 10-20 kpc in radius 
of NGC 5044. 
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Throughout this paper, we assume the Hubble constant to 
be Hq — lOO/i km s“^Mpc“^ and use the 90% conhdence 
level. At the redshift of 0.0090 (from the he NASA/IPAC 
Extragalactic Database; NED) for NGC 5044, 1' corresponds 
to 7.7h~^ kpc. We use 4.68 x 10“^ in number for the solar Ee 
abundance. 


2. Observations 


XMM-Newton observations of NGC 5044 were performed on 
2001 January 12. The RGS dispersion axis was oriented at a 
position angle of 114 degrees (North to East). We obtained a 
useful RGS exposure time of 22 ksec. Detailed descriptions of 
the XMM-Newton satellite and the RGS instrument are found in 


Jansen et al. (2001) and den Herder et al. (2001), respectively. 
Buote et al. (2002) analyzed the EPIC data of this observation 
and submitted a paper on the spatial properties of the X-ray 
emitting gas. 


3. Analysis and Results 

3.1. Data extraction and fitting method 


We analyzed the R GS data in a way similar to previous reports 
(e.g. Tamura et al. 2001b ), but with improved calibration of 
the instrumental responses and background. Eor basic data pro¬ 
cessing, we used the Science Analysis System (SAS ver.5.3.3). 
We extracted two sets of spectra from two different regions, 
|r(;| < 1' and 1' < |r(;| < 2', where w is the distance from the 
X-ray maximum in the cross-dispersion direction. Hereafter we 
refer to these spectra as W02 and W24, respectively. Because 
of the principle of the instrument, each spectrum contains emis¬ 
sion from a range of different projected positions. Based on the 
observed EPIC/MOS image in the RGS energy band, we es¬ 
timate this mixing, as given in Table [^. This estimate implies 
that more than 80% of the flux originates from r < 2' and 
1' < r < 4' for W02 and W24, respectively, where r is the 
projected radius from the center. Since the instrumental spa¬ 
tial resolution, ~ 15" in EWHM, is smaller than our extraction 
widths, the PSP effects insignificantly on this estimation. 

We present the spectra (sum of W02 and W24) in Pig. |^. 
Thus, we can resolve lines from H- and/or He-like ions of Si, 
Mg, and O, along with a series of Pe-L lines. The observed 
Pe-L lines are from ionization states ranging from Pe XVII to 
PeXXl, suggesting emission from their maximum formation 
temperatures ranging from 0.4 keV to 1.1 keV. 

The emission is extended on larger scales than the instru¬ 
mental spatial resolution resulting in degradation of the line 
spread function (ESP). We approximated this effect by con¬ 
volving the source surface brightness p rofile w ith the RGS re¬ 
sponse for a point source (Kaastra et al. 2002a). In this way, the 
line width approximately becomes 


[f + {0.12 k ■ e/rf] 


210.5 


( 1 ) 


, where I and 9 are intrinsic instrumental line width and source 
spatial width, respectively. Here we modeled the brightness 
based on the observed MOS image. Due to the mirror PSP, 
the MOS image is more extended than the actual brightness 


Table 1. The estimated flux mixing effect, based on the MOSl 
image**". 


2D-radius (') 

W02 

W24 

<1 

0.28 

0 

1-2 

0.16 

0.15 

2-4 

0.07 

0.10 

4-5 

< 0.03 

< 0.04 


** The numbers give fractions of fluxes contributed by individual 
two-dimensional annuli to the RGS data integration regions, rela¬ 
tive to the total flux within 5* in radius. 



10 15 30 35 30 



Fig. 1. The RGS spectrum extracted from the central 4' in full- 
width from the NGC 5044 group. The spectrum is corrected 
for redshift and effective area, but not for the background. Pour 
spectra from the two instruments and two spectral orders (m=- 
1 and -2) are combined. The bottom panel is a zoom of a part 
of the spectra. The positions of H- and He-like ions (top panel) 
and Pe-P transitions (bottom) are indicated. 


profile. To compensate for this for the W02 spectral fitting, 
we re-scaled the brightness distribution to match the observed 
RGS profiles of the O VIII and Pe-P lines, observed in the 14- 
20.5 A range. We found a scaling factor in wavelength width 
of / = 0.83, which is used for the spectral fitting below. This 
corresponds to 9 = (observed width with MOS) x/. Por W24, 
which has too poor statistics to determine the scaling factor, we 
simply used / = 1. 

We limited the analysis to the wavelength band of 8-28 A, 
where the estimated background is typically less than 50% of 
the source flux. The background spectrum and its variability 
were estimated from several blank-sky field observations (~ 
400 ks in total exposure). A 30% systematic error was assigned 
to the background counts (den Herder et al. 2002). The first 
order spectra from the RGSl and RGS2 were fitted with the 
same model simultaneously. 
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3.2. Isothermal model fitting 


To characterize the observed spectra and hnd the temperature 
and metallicity, we htted the spectra with a collisonal ionization 
equilibrium model (CIE; implemented in the SPEX package; 
Kaastra et al. 2002a) modihed by photoelectric absorption and 
cosmological redshift (IT model). The absorbing column den¬ 
sity was hxed to the Galactic value, 4.9 x 10^"* m“^, obtained 


from the Hi map (Dickey & Lockman |1990| using NASA’s 
W3nH tool). The abundance ratios relative to the solar value of 
He was hxed to unity (He/H=l), while those of other heavy ele¬ 
ments, except for O, were hxed to that of Ee (e.g., Si/Ee=l). The 
solar abundances are taken from Anders and Grevesse (1989). 


The free parameters are the emission measure (EM), the tem¬ 
perature, Ee/H, and O/Ee ratios. 

The results of htting with the IT model are shown in 
Table ^ and Eigures |[^. Although the hts are not acceptable, 
the IT model reproduces the global spectral structure fairly 
well, indicating that the emission is dominated by a single tem¬ 
perature component. The representative temperatures derived 
from the IT ht, 0.8 keV (W02) and 0.9 keV (W24), are con¬ 
sistent with the p revious ROSAT and A SCA measurements 
(David et al. 1994; Eukazawa et al. 1996). 


Signihcant residuals are found around 11-12 A, where 
lines from Ee XXII - Ee XXIV are expected, suggesting the pres¬ 
ence of emission hotter than the effective temperature. Another 
excess can be seen around 16 A in both the W02 and W24 spec¬ 
tra; the model underestimates the data. This is unlikely due to 
the incorrect modeling of the LSE, because the model repro¬ 
duces other lines, which should have similar brightness distri¬ 
bution and hence similar LSE. This is perhaps intrinsic to the 
emission model of lines (mainly Ee XVIII and O VIII Ly/j). In 
fact, the RGS spectrum of Capella also shows excess aroun d 
this band above the SPEX/CIE model (Audard et al. 2001 ). 
Another residual above the IT model can be seen around 22 A. 
Although this is close to the red-shifted position of O VII lines, 
we presume that this is due to residual instrumental hot pix¬ 
els/columns instead of emission from the source, because it is 
too narrow compared to other lines. There are some hot pix¬ 
els/columns in the RGS CCDs whose duty cycles are too low 
to detect statistically within normal integration time. 


3.3. Multi temperature model fitting 

Given the IT ht results, we attempt to obtain a better ht by 
introducing a combination of multiple CIE components with 
different temperatures. 

In the model htting, we hxed the temperature separation as 
Ti/Ti+i = F = constant, where Ti is the Ath temperature. 
A smaller F gives a better description of true EM structure, 
but at the same time it causes coupling between the compo¬ 
nents. Due to the nature of the Ee-L line emissivities, we can¬ 
not resolve the EM structure with resolution hner than about 
F = 1.25, regardless of the data quality (Kaastra et al. |2002b ). 
Given the present spectral quality, we use F ~ 1.5. We al¬ 
lowed each component to have a free normalization, EM{Ti), 
but hxed abundances among components. Therefore, additional 
free parameters are To and EM{Ti), where To is the tempera- 



Wavelength (A) 


Fig. 2. The RGS hrst order spectra extracted from a full-width 
of 2' (W02; hlled-circle for RGSl, cross for RGS2). Top, mid¬ 
dle, and bottom panels show the data compared with the best-ht 
IT model (histogram), the residuals (data-model/model) for the 
IT model, and those for the multi-T model, respectively. In the 
top panel, the background spectra are shown with dashed-line 
histograms. 
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Fig. 3. Same as the previous hgure, but for the W24 spectra. 


ture of the hottest component. We limited T^ to be hotter than 
0.1 keV since the RGS spectrum is not sensitive to much cooler 
emission. We used the best-ht Ee and O abundances obtained 
in the IT ht as initial parameters. Eurthermore, for W02 we 
started with Tq = 1.2 keV, which is 1.5 times the representa¬ 
tive temperature. Eor the W24 ht, we hxed Tq to the best-ht 
value of the W02 ht (1.07 keV) , because we are interested in 
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Table 2. The results of model fits to the RGS spectra. 


(1) 

(2) EM 

(3)T 

(4) Fe/H 

(5) O/Fe 

(6) AAh 

(7) X^lv 

model 

W02'‘ 

3.8 

0.78 

0.41 

0.54 

Of 

Aimi% 

IT 

W24^ 

2.0 

0.88 

0.38 

0.43 

Of 

284/183 

IT 

W02 

- 

1.07±0.03 

0.55±0.05 

0.51±0.06 

Of 

339/273 

multi-T 

W24 

- 

1.07f 

0.53±0.05 

0.42±0.12 

Of 

242/179 

multi-T 

W02 

- 

1.07f 

0.48±0.05 

0.57±0.06 

'j c + b.U 

'•'^-1.0 

353/279 

ICF-i-AAh 


(1) The spectral extraction region; W02 and W24 as described in the text. 

(2) The volume emission measure {h~^ 10^°m“®). 

(3) The temperature (keV). The highest temperature (To) in the case of a multi-T and ICF models. 

(6) Column density (10^® m“^). The covering factor of the excess absorption for the hot component is fixed to be 0.1. 

No errors are given for these fits, because the models are far from acceptable, 
f Fixed parameters. 


the difference in the EM shape between W02 and W24, in¬ 
stead of the value of Tq. 


As shown in Table ^ and Figures the fits improved sig¬ 
nificantly ( Ax^ =139 and 42 for W02 and W24, respectively). 
The W02 spectra can be described by the sum of two CIE com¬ 
ponents (To ~ 1.1 keV and Ti ^ 0.7 keV). A relatively small 
contribution from a 0.14 keV component is marginally sug¬ 
gested. No other component (0.2-0 .6 keV) is required as illus¬ 
trated in Fig. ^ where we compare the derived EM distribution 
with that predic ted fro m an isobaric cooling flow model (ICF; 
Canizares et al. |l988 ). The ICF model is the simplest descrip¬ 
tion of emission from a cooling gas without any heating; the 
differential EM (dEM / dT) is proportional to the inverse of 
the cooling function and the mass deposition rate. This compar¬ 
ison clearly indicates a lack of the emission from 0.2-0 .6 keV 
gas, relative to the emission from 0.7 keV gas, contrary to the 
prediction of the ICF model. 


A direct comparison of the observed spectrum with the ICF 
model shows that the 0.7 keV component nicely reproduces the 
fluxes and ratios of major resolved lines (Fe XVII, Fe XVIII, and 
O VIII) at least within a wavelength range of 14-20 A. Further, 
the 0.5 keV or 0.3 keV components of the ICF model clearly 
overpredict line fluxes of Fe XVII at 15.0 A and 17.0 A. Finally, 
in the 0.2 keV and 0.1 keV ICF model components, the line 
fluxes from O VIII and O VII become so high relative to the Fe- 
L lines that they are inconsistent with the data. 


In addition to O and Fe, we attempted to constrain the abun¬ 
dances of N, Ne, Mg, and Si. The H- and/or He-like lines 
of these elements are covered by the RGS. To derive indi¬ 
vidual abundances separately, we fixed all parameters except 
the relevant abundance and fitted within the limited spectral 
wavelength bands around the expected line from each ele¬ 
ment (Table ||). We used the best-fit multi-temperature model. 
Table shows the result for W02. Among these abundances, 
those of N and Si could have systematic errors larger than the 
statistical errors shown in the Table. In the case of N, emis¬ 
sion from gas cooler (<0.1 keV) than the RGS sensitivity may 
contribute to the line emission. The calibration of the RGS re¬ 
sponse around the Si line (^ 7 A) is relatively less certain; the 
Si abundance may be subject to systematic errors up to 20%, 
which is similar to the statistical ones. Therefore we conserva¬ 



Fig. 4. The constraint on the differential emission measure in 
terms of The filled circles and triangles indicate the best- 
fit values for the W02 and W24 spectra, respectively, while the 
open circles indicate the estimated values for the emission from 
the galaxy core (r < 1') only of the 0.7 keV and 1.1 keV com¬ 
ponents. The dashed and dotted lines indicate upper limits on 
values for W02 and W24, respectively. The model predicted 
from the isobaric cooling flow with a mass deposition rate of 
10 h~‘^ M 0 yr“^ is shown as the solid line. 


tively conclude that these elemental abundances relative to Fe 
are consistent with the solar ratio within 50% uncertainty. 

3.4. Absorption by cold material 

One possibility to hide the cool emission is absorption by blobs 
of cold material. If the cold material is located close to or as¬ 
sociated to the X-ray emitting cool component, the expected 
X-ray emission will be selectively absorbed. This cold material 
should also absorb a part of the hot component emission. In 
the ICF model, the relative volume occupied by a temperature 
component is proportional to dT / [A(T)n^] oc T^dT/A(T), 
where A(T) and n are the emissivity and the gas density. This 
implies that the volume occupied by the 0.4 keV component is 
roughly 10% of the volume of the 1 keV component. Therefore, 
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Table 3. The abundances of lighter elements relative to iron, in 
solar units'^. 



Abundance 

band^’lA) 

N/Fe 

0.62 ±0.5 

24-25 

Ne/Fe 

0.55 ±0.3 

11-14 

Mg/Fe 

1.26 ±0.2 

8-10 

Si/Fe 

1.50 ±0.3 

6-7.5 


Derived from the W02 spectra, using the best-fit multi¬ 
temperature model (Fe/H=0.55). 

^ The wavelength band used in the analysis. 


the covering fraction for the hot component to be absorbed 
is expected to be also ^ 10%. To examine this possibility, 
we modeled the RGS spectrum by the ICF model and differ¬ 
ential absorption. Assuming that 100% and 10% of the cool 
(0.1 < kBT/keV < 0.6) and hot (0.7 < fcsT/keV < 1.1) 
emission are absorbed by the same column density, we have 
confirmed that a column density higher than 6.5 x 10^® m~^ is 
required to describe the observed spectrum (Table Note that 
the best-fit Fe and O abundances did not change significantly 
in this modeling. 

If there is cold material with A^h — 10^® m“^, more than 
99 % of the soft X-rays (< 0.7 keV) must be absorbed, re¬ 
sulting in emission holes for external observers. To examine 
whether there exist such brightness depressions, we have ana¬ 
lyzed archival Chandra data of NGC 5044. As shown in Fig. ||, 
there is no such effect near the center, nor drastic inhomogene¬ 
ity of the emission hardness, down to a spatial scale of ~ 10" 
(~ 2 kpc). Therefore, we conclude that this scenario is very 
unlikely. 


4. Discussion 

4 . 1. Summary of the results 

We have successfully resolved a number of emission lines from 
the central regions of the NGC 5044 group. As a result, we have 
found that the spectra extracted with a full-width of 2' (W02) 
can be described by a combination of a 0.7 keV and 1.1 keV 
CIE component. Any cool component with temperature in the 
range of 0.2-0.5 keV and 0.1-0.2 keV is less than a few and ten 
percent of the overall emission measure (EM), respectively. 

4.2. Spatial temperature distribution 

The W02 spectra contain emission from different projected 
radii. We can approximately subtract the contribution from the 
outer region by assuming that the spectral shape is azimuthally 
uniform and using the surface brightness distribution observed 
with the MOS. The estimate in Table implies that similar 
fluxes from a region of 1' < r < 4', where r is the pro¬ 
jected radius, are contaminated into the W02 and W24 spectra. 
Therefore, by subtracting the W24 model from that of W02 we 
can estimate the spectral nature of the galaxy core (r < 1'). As 
shown by open circles in Fig. ^ the projected core (and hence 
the three-dimensional core) is dominated by the 0.7 keV com¬ 


NGC 5044 (0.3 - 1 keV ) 



SKY X (ai-csec) 



SKY X (arcsec) 


Fig. 5. Chandra ACIS-S image of the center of NGC 5044. Top 
is the 0.3-1.0 keV brightness map smoothed with a Gaussian 
of cr = 1".5, and bottom is a hardness ratio map (1-2 keV / 
0.3-1 keV) of the same region. The scales are arbitrary in both 
panels. 


ponent and the EM of the 1.1 keV component is less than 20% 
of that of the 0.7 keV emission. 


This RGS result is consistent with the EPIC results by 
Buote et al. (2002). They deprojected the observed spectra into 
three dimensional space and found that within and beyond 1' 
the temperature is 0.75 ± 0.1 keV and 1.2 ± 0.1 keV, respec¬ 
tively. 


Because the separation between the two temperatures de¬ 
rived from the RGS (0.7 keV and 1.1 keV) is small, it is diffi¬ 
cult to distinguish between (1) emission from two temperature 
components and (2) a continuous distribution in temperature 
between the two temperatures. The similarity in the obtained 
temperature values between the RGS and EPIC, which have 
different energy responses, suggests that model (1) is more 
likely than model (2). 
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4.3. Possibilities of resonance line scattering 

Resonance line scattering could affect line strengths from the 
core regions of X-ray bright clus ters a nd elliptical galaxies 
(Gil’fanov et al. 1987 ; Shigeyama 1998 ). For the temperature 
of the core of NGC 5044, ~ 0.7 keV, the highest optical depth is 
expected for F e XVII at 15.0 A (oscillator strength / = 2.7). In 
fact, Xu et al. (2001) detected resonant scattering effects in the 
15.0 A line, from the RGS observations of the core of the ellip¬ 
tical galaxy NGC 4636. This galaxy has a temperature structure 
similar to that of NGC 5044. 

In the case of NGC 5044, the optical depth of this line to¬ 
wards the galaxy center is estimated to be about 4 when there 
is no line broading due to bulk motion of the gas. Then, the 
15.0 A line photons would be scattered out of the core region, 
resulting in a depression in W02 and an enhancement in W24. 
However, there is no evidence for this in the observed spec¬ 
tra (Figures This is probably because we have integrated 
all emission within a full-width of 2', which is wider than the 
emission core of this object, ~ 30" in radius. Thus, we estimate 
the resonance scattering effects to be small. 


4.4. Origin of the lack of cool component 

According to our results and the above argument, we presume 
that there is a cut-off in temperature, T^ut, in the true EM dis¬ 
tribution. In cores of medium-size to rich clusters, a similar 
lack of cool emission has been reported based on the high res¬ 


for elliptical galaxies (e.g. Matsumoto et al. 1997; Matsushita 
2001): fj,mpa^/kT = 0.5 — 1.0, where fi ~ 0.6 and rup are 


the mean molecular weight and proton mass. Furthermore, by 
fitting the ASCA spectra of ^ 10 cluster centers by a two tem¬ 
perature model, Ikebe (2001) found a tight correlation between 
the two temperatures, which resembles the present correlation 
between Tcut and the ICM temperature (Fig. ^c). It is unclear, 
however, whether these correlations are intrinsic or due to the 
correlation of Tcut with the ISM temperature combined with 
the correlation of the ISM temperature with these parameters. 

Based on the correlation between Tcut and the scale of 
the galaxy, at least a substantial part of the heating energy 
may originate from the galactic gravitational potential. Possible 
energy-transfer mechanisms include an adiab atic co mpression 
along with the gas inflow (e.g. David et al. 1994) and stel¬ 
lar motion through mass-loss or amplification of interstellar 
magnetic fields (e.g. Ikebe et al. 1999). These heating effects 
should be taken into account in relatively smaller systems like 
NGC 5044, where the gravitational energy is comparable to the 
thermal energy of the gas. For the larger systems, however, the 
gravitational energy may be not enough to balance the large 
X-ray luminosity. Alternatively, this correlation may originate 
from a direct interaction between the gas and the dark matter 
(Qin& WupMib. 


olution RGS spectra (Kaastra et al. 2001; Peterson et al. 2001 


Tamura et al. 2001a, 2001b; Sakelliou et al. 2002, Xu et al. 


2001), as well as from A5CA (e.g. Makishima et al. 2001) and 
Chandra (e.g. David et al. 2001). Among these, our result along 


The correlation between Tcut and the ICM temperature and 
density also suggests that the surrounding environment is also 
an important parameter to determine Tcut- We further specu¬ 
late that a thermal coupling between the central region and 
the ICM has been taking place. For example, a certain mag¬ 
netic field configuration may transfer a large amount of energy 
from the ICM to the galactic cool gas as the member galaxies 


with that of NGC 4636 by Xu et al.(|2001|) is unique because move through the ICM (Makishima et al [200 1| and references 


1) the quality of the observed spectra is among the best, 2) 
the observed Tcut value is significantly lower than for other 
sources studied, and 3) the galaxy is surrounded by a relatively 
cooler (and lower pressure) intracluster medium (ICM) than 
in other cases. A comparison with the ICF prediction suggests 
that some mechanism has heated or mixed the gas with a to¬ 
tal energy (~ 10^® W in 10 Gyr) and efficiency high enough 
to counterbalance the observed X-ray radiative cooling in this 
system. Several ideas for the heatin g and mixing have b een dis¬ 
cussed)^. Makishima et al. 2001 ; Peterson et al. POOlt Fabian 


therein). Alternatively, the jet from active galactic nuclei at the 
galaxy center and associated bubbles may stir the cool g alacti c 
gas into surrounding hot ICM gas (e.g. Churazov et al. ^001 ). 
Note, however, that among six systems shown in Fig. ^ only 
M87 exhibits a well dehned radio jet or bub ble. NGC 5044 has 
only a compact radio core (Slee et al. 1994 ). One important test 
for this model is to check if this large scale mixing preserves the 
observed metallicity concentration around the cD and other gi¬ 


et al. 2001 


Brighenti & Mathews 2002; Bohringeret al. 2002). 


ant galaxies (e. g. Fuji ta & Kodama |1995| ; Fukazawa et al. |1998 
Allen & Fabian 


To further understand the possible heating or mixing in 
these systems, we compare the values of Tcut observed from 
several giant elliptical galaxies, against other observables as 
shown in Fig.^. These galaxies are in the outskirts (NGC 4636) 
or center (M 87) of the Virgo cluster, in the center of a group 
(NGC 5044), and in the centers of poor (Sersic 159-03) to rich 
(A 496, A 1795 and A 1835) clusters. The value of Tcut was de¬ 
termined for all cases from the RGS spectra, but not uniformly 
analyzed and therefore include a systematic uncertainty up to 
30%. As is clear from the figure, the values of Tcut of these 
galaxies appear to be correlated to their optical luminosity and 
the stellar velocity dispersion (cr), as well as to the properties of 
the surrounding ICM. Note that a correlation between the the 
emission-weighted gas temperature and a has been reported 


To investigate each model and understand the physics 
around central galaxies, systematic observations of galaxies 
on different scales and in different environments along with 
theoretical studies are crucial. For such attempts, we refer to 


Kaastra et al. (2002b) and Peterson et al. (2002[). 


4.5. Origins of metals in galaxies and intra-cluster 
medium 

Our accurate measurement of the thermal structure has in turn 
provided the determination of the Fe and O abundances, more 
robustly than any previous measurements. Using the EPIC data 


of NGC 5044, Buote et al. (2002) reported that the Fe abun¬ 


dance (R < 2') is 1 ± 0.2 relative to the meteoritic solar abun- 
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dance[|. This value corresponds to ~ 0.7 ± 0.15 solar in our 
definition, and is consistent with our measurement, 0.55±0.05, 
within the errors. 

The Fe abundance of NG C 5044 we measured and that of 
NGC 4636 (Xu et al. 2001), 0.5-1.0 times solar, are within 
the range of ASCA measurements in 27 ellipticals (Matsushita 
et al. 2000| ). Therefore we refer the reader to the discussion 
on the Fe abundances presented by Matsushita et al. (2000). 
Most importantly, we have confirmed that the Fe abundance 
in the X-ray emitting gas is already comparable to that of the 
stellar component and consequently additional enrichment by 
supernovae is not substantially required, when no dilution by 
the in-flowing ICM is assumed. 

Contrary to Fe, the determination of the O abundance in 
NGC 5044 is reported here for the first time. The O/Fe ratio 
of NGC 5044, about half of the solar ratio, is similar to that of 
NGC 4636 and other galaxies at the cluster core (e.g., Tamura 
et al. 2001b; Kaastra et al. 2001), indicating a common origin 
of the metals in these systems. On the other hand, the O/Fe ra¬ 
tios in these ellipticals are about two times smaller than those 
in the intracluster medium (e.g. see the previous references) 
and the solar ratio. This difference suggests the presence of at 
least two different origins for metals in the ellipticals and the 
ICM, regardless of any theoretical models. Variati ons in the 
Si/Fe ratio within a cluster (e.g. Fi noguen ov et al. 2000 ) and 
among clusters (e.g. Fukazawa et al. 1998 ) also support a mul- 
tipl e origin of the metals in the ICM. For example, Tamura et 
al. ( 2001b ) explained the change in the O/Fe ratio between the 
central galaxy and the ICM in A496, in terms of an increase of 
the relative contribution from type la supernovae to the metal 
production in the galaxy. The present result supports this sim¬ 
ple idea. 

To further understand the metal production in clusters, sys¬ 
tematic measurements of variations in the O/Fe ratio within a 
cluster and among clusters are important. Note that the O/Fe 
ratio is more sensitive to their origins than the Si/Fe ratio. 
For such a study the RGS has not enough sensitivity for spa¬ 
tially extended O lines. Future experiments such as X-ray mi¬ 
crocalorimeter (XRS) on board Astro-E II would provide more 
accurate measurements. In addition to the O, Si, and Fe abun¬ 
dances, the C and N abundances are also valuable. This is be¬ 
cause C and N are believed to be produced largely by stellar 
mass loss instead of supernovae. To measure those abundances 
a higher sensitivity for lower energy X-ray lines (0.1-0.5 keV) 
is needed. 
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